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1 Introduction 



Quantum Chromodynamics (QCD) is regarded as the correct theory for 
strong interactions. Investigation of bound states of quarks (and/or glu- 
ons) is one of the most effective methods to study this dynamics among 
these constituents. Since the task of calculating hadron structure from QCD 
itself is very difficult (as can be seen from various lattice QCD methods), one 
can on the other hand rely on specific models to gain some understanding 
of it at low energies, and this study can most effectively be accomplished 
by applying a particular framework to a diverse range of phenomena. Me- 
son decays provide an important opportunity for exploring the structure of 
these simplest bound states in QCD and for studying the non-perturbative 
(long distance) aspect of the strong interactions. Besides vector mesons (see, 
e.g., P] and refs. therein), pseudo-scalar mesons have for a long time been a 
major focus of attention to understand the inner structure of hadrons from 
non-perturbative QCD. A number of such studies [2J El IU El El Ej dealing 
with decays of pseudo-scalar mesons at quark level of compositeness have 
been carried out recently. In this paper, we study leptonic decays of pseudo- 
scalar mesons such as K, D, D$ and B, which proceed through the coupling 
of the quark-antiquark loop to the axial vector current as shown in Figure 1. 

A realistic description of pseudo-scalar mesons at the quark level of com- 
positeness would be an important element in our understanding of hadron 
dynamics and reaction processes at scales where QCD degrees of freedom 
are relevant. Such studies offer a direct probe of hadron structure and help 
in revealing some aspects of the underlying quark-gluon dynamics. The rel- 
ativistic description for analyzing mesons as composite objects is provided 
by the framework of Bethe-Salpeter Equation (BSE) in this paper. We em- 
ploy the QCD oriented BSE under Covariant Instantaneous Ansatz (CIA) 
[8]. CIA is a Lorentz-covariant generalization of Instantaneous Approxima- 
tion. For qq system, CIA formulation [8] ensures an exact interconnection 
between 3D and 4D forms of the BSE. The 3D form of BSE serves for mak- 
ing contact with the mass spectra, whereas the 4D form provides the Hqq 
vertex function T(q) for the evaluation of various transition amplitudes. A 
BSE framework under Instantaneous Approximation formulation similar to 
the CIA formulation was also earlier suggested by the Bonn group [9]. 

We had earlier employed the framework of BSE under CIA for calcula- 
tion of decay constants [HI [TO] of heavy-light pseudo-scalar mesons and F n for 
7T° — > 27 process. We also evaluated the leptonic decays of vector mesons, 
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such as p, u, 4> [H] in this framework. However, one of the simplified assump- 
tions in all these calculations was that the Hqq vertex was restricted to have 
a single Dirac structure (e.g., 75 for pseudo-scalar mesons, 7 • e for vector 
mesons, etc.). However, recent studies [121 13] have revealed that vari- 
ous mesons have many different covariant structures in their wave functions 
whose inclusion was also found necessary to obtain quantitatively accurate 
observables [12J and it was further noticed that all Dirac covariants do not 
contribute equally and only some of them are relevant for calculation of me- 
son mass spectrum and decay constants. Such a copious Dirac structure 
of the BS wave function in fact was already indicated by Llewellyn Smith 
Hence it is necessary to introduce various Dirac structures into the Hqq 
vertex for different kinds of mesons. In the recent work pQ, we developed a 
power counting rule for incorporating various Dirac covariants in the struc- 
ture of vertex function, order by order in powers of inverse of meson mass, 
and calculated the leptonic decays of equal mass vector mesons such as p, u, 
<f), taking into account the leading order covariants since they are expected to 
contribute maximum to observables according to our scheme. On the line of 
that work, in this paper we first discuss the power counting rule for choosing 
various Dirac covariants from their complete set (see, e.g., [51 [121 131 Ej) for 
pseudo-scalar mesons in Section 2. In section 3 we calculate leptonic decay 
constants of them employing the wave function developed in section 2. We 
then conclude with discussions in Section 4. 

2 Structure of generalized vertex function T(q) 
for pseudo-scalar mesons in BSE under CIA 

For introducing the variables and for the convenience to discuss the gen- 
eralized hadron-quark vertex in BS wave function under CIA, we give the 
similar outline for CIA as in Ref. [H [11]. We start with a 4D BSE for scalar 
qq system with an effective kernel K and 4D wave function $(P, q): 



2(27r) 4 A 1 A 2 $(P,g) = J d 4 q' K(q, q')$(P, q') 



(1) 



where A 12 , the inverse propagators of two scalar quarks, are: 



A lj2 = m\ 2 +pl 2 . 



(2) 
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Here m 12 are (effective) constituent masses of quarks. The 4-momenta of 
the quark and anti-quark, p 12 , are related to the internal 4-momentum g M 
and total momentum P of hadron of mass M as 



Pi,2fM = rhx^Pfi ± <? M , (3) 

where = [1 ± (m| — m2)/M 2 ]/2 are the Wight man- Gar ding (WG) defi- 
nitions [10] of masses of individual quarks. 

The CIA Ansatz on the BS kernel K in Eq. (JTJ) is, 

K(q,q')=K(q,q% (4) 

where 

g • P 

% = %- -p2~Pn (5) 

is observed to be orthogonal to the total 4-momentum P ( i.e., q.P = 0), 
irrespective of whether the individual quarks are on-shell or off-shell. A 
similar form of the BS kernel was also suggested in ref. [9]. The longitudinal 
component of q^, 

Ma = M q -^, (6) 

does not appear in the form K(q, q') of the kernel. For reducing Eq. (JTJ) to 
the 3D form, one can define a 3D wave function <p(q) as 

+oo 

0(g) = f Mda$(P,q). (7) 

— oo 

Following usual steps outlined in [TJ [8j [TTJ , we get the Hqq vertex function 
T(q) under CIA for the case of scalar quarks: 

— oo 

where D(q) is a 3D denominator function whose value can be easily worked 
out by contour integration by noting the positions of the poles in the complex 
a— plane [TJ [8j [TTJ. By this process, an exact interconnection between 3D 
wave function <fi(q) and the 4D wave function $(P, q) and hence that between 
3D and 4D BSE is thus brought out, where the 3D form serves for making 
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contact with the mass spectrum of hadrons, whereas the 4D form provides the 
Hqq vertex function T(q) which satisfies a 4D BSE with a natural off-shell 
extension over the entire 4D space (due to the positive definiteness of the 
quantity q 2 = q 2 — ^ q 'D throughout the entire 4D space) and thus provides a 
fully Lorentz-covariant basis for evaluation of various transition amplitudes 
through various quark loop diagrams (see Figure 1). 

To apply the above simplified discussions to the case of fermionic quarks 
constituting a particular meson we proceed in the same manner as p]: The 
scalar propagators in the above equations are replaced by the proper 
fermionic propagators Sp. Then, on observing the vertex T(q) now is a 4 x 4 
matrix in the spinor space, we should incorporate its relevant Dirac struc- 
tures, for which we take guidance from some recent studies [3 [T2J, [13j|, as 
well as Llewelyn-Smith's classic paper [TJ], which have revealed that vari- 
ous mesons have many different covariant structures in their wave functions 
whose inclusion was found necessary to obtain quantitatively accurate ob- 
servables. It was also noticed recently [12] that all Dirac covariants do not 
contribute equally and only some covariants are considered to be relevant for 
calculation of mass spectrum and decay constants. Ref. [12] has also calcu- 
lated masses and decay constants of vector mesons for various subsets of co- 
variants from their complete set. Towards this end, we make use of the power 
counting rule developed in [T] for incorporating various Dirac covariants in 
the structure of Hqq vertex function for a particular meson, order-by-order 
in powers of inverse of the meson mass M, so as to systematically choose 
among various covariants from their complete set and write wave functions 
for various mesons. 

As far as a pseudo-scalar meson is concerned, its Hqq vertex function, 
which has a certain dimensionality of mass, can be expressed as a linear 
combination of four Dirac covariants Tf(i = 0, 3) [121 DSl EH]) each mul- 
tiplying a Lorentz scalar amplitude Fi(q 2 , q ■ P, P 2 ). The choice of the Dirac 
covariants is not unique as can also be seen from the choice of covariants 
used in Ref. [T21CI3]. For adapting this decomposition to write the structure 
of vertex function T(q), we re-express the Hqq vertex function by making 
the amplitudes Fi(q 2 , q ■ P, P 2 ) dimensionless, weighing each Dirac covariant 
with an appropriate power of M. Thus each term in the expansion of T(q) 
is associated with a certain power of M. In detail, we can express Tp as a 
polynomial in various powers of 1/M: 

T P = n p J-N P D(q)cl>(q), (9) 
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with 



B B B- 

tip = 75.00 -275(7 - p )~^[ -^5(7- -75(7- Pi -g-7- 11 ' P ">J^2- ( 10 ) 

Here Bi(i = 0, ...,3) are four dimensionless and constant coefficients (which 
are taken to be constant on lines of PQ, only to consider the leading powers of 
1/M) to be determined. Now since we use constituent quark masses where 
the quark mass m is approximately half of the hadron mass M, we can use 
the ansatz 

q«P~M (11) 

in the rest frame of the hadron (among all the pseudo-scalar mesons, pion 
enjoys the special status in view of its unusually small mass (< Aq C d) and its 
case should be considered separately. See the discussions in Section 4). Then 
each of the four terms in Eq. receives suppression by different powers of 
1/M. Thus we can arrange these terms as an expansion in powers of O(j^). 
We can see in the expansion of Qp that the structures associated with the 
coefficients B ,Bi have magnitudes O(^j) and are of leading order, while 
those with B2, -B3 are 0(-^j-)and are next-to-leading order. This naive power 
counting rule suggests that the maximum contribution to the calculation 
of any pseudo-scalar meson observable should come from the leading order 
Dirac structures 75 and ij^.P)-^ associated with the constant coefficients 
B and Bi, respectively. As a first application of this to the pseudo-scalar 
meson sector and on lines similar to [lj for vector meson case, we take the 
form of the Hqq vertex function incorporating the leading order terms in 
expansion ([9]) and ignoring O(tjt) terms for the moment, i.e., 

r(g) = [ 75 5o - 2757 • p^)^N P D{q)4>{q) ■ ( 12 ) 

As has been stated in p], the restriction from charge Parity on the wave 
function of the eigenstate should be respected 0. 

From the above analysis of the structure of Hqq vertex function in Eq. (ll2l) 
we notice that, at leading order, the structure of 3D wave function 0(g) as well 
as the form of the 3D BSE are left untouched and have the same form as in our 



1 To get the complete set of the Dirac structure for a certain kind of mesons, the restric- 
tion by the (space) Parity have been employed; and it is easy to see that the requirements 
of the space Parity and the charge Parity are the same for the vertex as well as the full 
wave function (see also [14]). 
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previous works, which justifies the usage of the same form of the input kernel 
we used earlier. Now we briefly mention some features of the BS formulation 
employed. The structure of BSE is characterized by a single effective kernel 
arising out of a four-fermion Lagrangian in the Nambu-Jonalasino [T5l [T6] 
sense. The formalism is fully consistent with Nambu- Jona-Lasino [16] picture 
of chiral symmetry breaking but is additionally Lorentz-invariant because of 
the unique properties of the quantity q 2 , which is positive definite throughout 
the entire 4D space. The input kernel K(q, q') in BSE is taken as one-gluon- 
exchange like as regards color (|A.|Aa) and spin ("f^^j 2 ^) dependence. The 
scalar function V(q — q') is a sum of one-gluon exchange Voce and a confining 

term Vb™/. m MH5]. Thus 

if(ff,^ = ^ 1) i^ 1 V«y( ff -g'); 

[q — q < Y 4 yy J uJq 

(13) 

M> = Max(M, mi + m 2 ). (14) 

The ansatz employed for the spring constant uj 2 ^ in the above equation is 
[1,11,15] 

u 2 q7j = 4:m 1 m 2 M > u; 2 a s (Q 2 ). (15) 

Here the proportionality of u 2 -^ on as(Q 2 ) is needed to provide a more 
direct QCD motivation to confinement. This assumption further facilitates 
a flavour variation in u 2 ^. And in Eq.(13) and Eq.(14) is postulated as a 
universal spring constant which is common to all flavours. 

In the expression for V(q — q'), as far as the integrand of confining term 
Vconf. is concerned, the constant term Co/cUq is designed to take into account 
the correct zero point energies, while term {clq « 1) simulates an effect 
of an almost linear confinement for heavy quark sectors (large m\,m2 ), re- 
taining the harmonic form for light quark sectors (small m 1; 7712) [HIS] , as is 



12tt 



33-2/ 



A 2 



7 



believed to be true for QCD. Hence the term r 2 (1 + 4a mim 2 M'^r' 2 )~2 in the 
above expression is responsible for effecting a smooth transition from har- 
monic (qq) to linear (QQ) confinement. The values of basic input parameters 
of the model are a = .028, C = .29, u = .158 (GeV) and A = .20GeV 
[1,10,11,15] which have been calibrated to fit the qq hadron mass spectrum 
obtained by solving the 3D BSE [15]. 

Now comes to the problem of the 3D BS wave function. The ground state 
wave function 0(g) satisfies the 3D BSE [lj on the surface P.q = 0, which 
is appropriate for making contact with 0(3)-like mass spectrum [15]. Its 
fuller structure (described in Ref. [15]) is reducible to that of a 3D harmonic 
oscillator with coefficients dependent on the hadron mass M and the total 
quantum number N. The ground state wave function 0(g)deducible from this 
equation thus has a Gaussian structure [U [HI [11] and is expressible as: 

0(g) ^e- ? / 2/32 . (16) 

In the structure of 0(g) in Eq. ffTBT) . the parameter j3 is the inverse range 
parameter which incorporates the content of BS dynamics and is dependent 
on the input kernel K(q, q'). The structure of j3 is given in Section 3. 



3 Decays constants fp of Pseudo-scalar Mesons 

Decay constants fp can be evaluated through the loop diagram shown in 
Figure 1 which gives the coupling of the two-quark loop to the axial vector 
current and can be evaluated as: 

f P P, =< 0\Qi lfll5 Q\P(P) >, (17) 

which can in turn be expressed loop integral: 

/pP, = V3|ATr[*p(P,?)%75] • (18) 
Bethe-Salpeter wave function \&(P, q) for a P-meson is expressed as 

(P, q) = S F ( Pl )r(q)S F (-p 2 ), with 

at \ .(m1-i7-.P1) Q ( \ .( m 2 + il-P2) nn , 
S F {Pi) = -i t ; S F {-p 2 ) = r • (19) 

In the following calculation, we only take the leading order terms in the 
structure of hadron-quark vertex function T(g) as in Eq. (fT2"|) . Sp are the 
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fermionic propagators for the two constituent quarks of the hadron and the 
non-perturbative aspects enter through the T(q). Using *$?(P, q) from Eq. (IT9l 
and the structure of Hqq vertex T(q) from Eq. (1121) . evaluating trace over 7- 
matrices and multiplying both sides of Eq. (1181) by P^j (— M 2 ), we can express 
fp as: 

f P = VZN P J d 3 qD(q)<f>(q)I , 

+00 9 
r Mda ,„„ . ,„ bm z . „. , , . 

' = /^A^ {2B » [m - (1 -MT ) + 2W| (2 °» 

— OO 

+ 5-[(-M 4 - 4 mi m 2 M 2 + m\ 2 5m 2 ) + 4M 2 g 2 - 4M 2 m 12 6ma + 4MV]}, 

where according to Eq.flS]) and Eq.([SJ, we had expressed scalar products prP2, 
Pi ■ Pand P2 ■ P as 

p 1 .p 2 = -M 2 (m 1 + ( x)(m 2 - ( T)-g 2 , (21) 

Px-P = -M 2 (mx + a) 
p 2 .p = -M 2 (m 2 -a), 

with m 12 = mi + m 2 and <5m = m 2 — mi. Here we have employed unequal 
mass kinematics when the hadron constituents have different masses. We 
see that on the right hand side of the expression for fp, each of the expres- 
sions multiplying the constant parameters Bo and B\ consist of two parts, 
of which only the second part explicitly involves the off-shell parameter a 
(that is terms involving a multiplying B , and terms involving both a and 
a 2 multiplying B\)). It is seen that the off-shell contribution which vanishes 
for mi = m 2 in case of fp calculation in CIA [10] using only the covariant 
75, now no longer vanishes for mi = m 2 in the above calculation, due to the 
term 4M 4 a 2 multiplying B\ when another leading order covariant ij^j-P/M 
is also incorporated in Hqq vertex function. This possibly implies that the 
off-shell part of fp does not arise from unequal mass kinematics alone, which 
is in complete contrast to the earlier CIA calculation of fp [10] employing 
only 75. This serves as a clear pointer to the fact that in this BS-CIA model, 
the Dirac covariants other than 75 might also be important for the study of 
processes involving large q 2 (off-shell), which is also suggested in [T7] . 
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Carrying out integration over da by noting the pole positions in the com- 
plex cx-plane 

Ai = 0=>ot = ±jjj[-m 1 T ie, u\ = m 2 + q 2 

A 2 = = +m 2 ±ie,u 2 = m 2 + g 2 , 1 ' 

we can express fp as 

// SttI/^ 1 

rf 3 gD(g)0(g)(2S o [m 12 (l - ^f)^py + ^mR,} 



+ — - [(— M - Am 1 m 2 M 2 + m 2 2 5m : 



M 3LV 1 z 12 

+ ^Ty(-4M 2 m 12 5m./? 1 +4M 4 J R 2 )]), (23) 

where the relationship between the functions A)(<z) and D(q) (see Ref. [H] 
for details) is D(q) = D (q) = {uj\ + uj 2 ) 2 — M 2 . The results of 

a contour integration is: 



Mda _ APj-UH + u 2 ) + (mj - m 2 )^ + u 2 ) 
>J ' 2mA! A/ 4M 2 u;iu; 2 (M 2 -(cu 1 + cc> 2 ) 2 ) ' 1 j 



"00 - ^ r 



27riAiA 2 

M 4 — m 2 l2 5m 2 + 4M 2 cj 1 c<j 2 )(c<j 1 + cu 2 ) + 2M 2 m 12 5m(w 2 — 
8MVc; 2 (M 2 - (cji + cj 2 ) 2 ) 



25) 



The structure of the parameter (3 in 0(g) appearing in Eq. fT23l is taken 
from Ref. [H DEI DH E51 HE! • (for details see Ref. pQ). Thus we take 



P 2 = (2m 1 m 2 Ma; g V 7 2 ) 1 /^ 7 2 = l - (26) 



Here is expressed as in Eq. f TToT) . To calculate the normalization f actor N P , 
we use the current conservation condition [8], 



2zP/i = (2vr) 4 j rf 4 gTr^(P,g)(^-^ 1 (p 1 ))^(P,g)^ 1 (-p 2 )] + (1 & 2), 



(27) 
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where the momentum of constituent quarks can be expressed as in Eq.flS]). 
Taking the derivatives of inverse of propagators of constituent quarks with 
respect to the total 4-momentum P M , evaluating trace over the 7-matrices 
and following usual steps outlined in Ref . [H [18] , then carrying out the da 
integral by noting the pole positions in complex cr-plane, we can express the 
normalizer as, 

Np 2 = -{2n) 2 i J d 3 qD 2 {q)<j) 2 {q)S ; with 

S = ^^l[(-s m 3 m 2 12 - M 4 m 12 + 2M 2 m 2 12 5m)I 1 

+ (2M 4 m 12 + 2M 4 5m + AM 2 m 12 5m 2 )I 2 - (AM 4 5m)I 3 ] 

- H [D^Exh + (2M 2 E 1 - AM 2 5mm 12 D 1 )I 2 

+ (-8M 4 5mm 12 + AM 4 D 1 )h 

+ (8M 6 )/ 4 - ^[D x Ekh + (2M 2 E 2 - AM 2 5mm l2 D l )I 2 

+ (-8M 4 5mm 12 + 4M 4 D 1 )/ 3 + (8M 6 )/ 4 ] 
ABoBi . „ 1 



+ —^\-2M 2 m 2 - 



M 3 1 D(qy 

r>2 1 

1 r /i 71 <r4 d 1 / oA/r2_ r„, 1 o/\/r4\ 



■ [AM R x + (-2M 2 mi 2 5m + 2M 



M 41 v y D(g) J 

||[-2(M 2 m 12 5m + M 4 )-^ + 4M 4 i? 1 ] 

q 2 (^l[AM 2 D 1 I 1 + 8M 4 / 2 ] - ^[-4M 2 J D 1 J 1 - 8M 4 / 2 ]).(28) 



M 4 l M 4 
In the above, the quantities D\,Ei,E 2 are: D\ = —ni\ 2 bm + M 2 ; E\ = 

8m 2 m\ 2 — AM 2 mim 2 — M 4 ; P 2 = — 3M 2 m 1 m 2 , while the integrals I\,I 2 , 13 
and I4 over the off-shell parameter da are: 



r _ / 



Mda 

2ni 



A, 



2u;f — M 2 cj 2 + 5tu 2 u; 2 + 4co>iu;f + cjf 
Aufu 2 (M 2 - (wi + cu 2 ) 2 ) 2 



(29) 



7 Mrfa 

J 2 = / 



A 2 A 



2 
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, -M A u 2 + (ml-ml)(u 1 + u 2 ) 2 (2^i + ^ 2 ) 

ZTlt 



8M 2 ufu 2 (M 2 - (ui + u 2/ 



2\2 



M 2 [Qm'f + 9Lujuj 2 + 4tJioj + uo 2 (-m 2 + m 2 + up] v 
+ 8M 2 c^ 2 (M 2 -( Wl +cu 2 ) 2 ) 2 ( j 



+ 00 > , T 



/3 = ' AfA^ 



— oo 

2M(M-m!)(M 2 + 



/ 2(M 2 -m 2 + m 2 + 2Mcj 2 ) 2 
Vcu 2 (M - U) X + u 2 ) 2 {M + u x + cu 2 ) 
- m!)(M 2 + m 2 - m 2 - 2Mcu 1 ) 2 



^((M-uj^-uj 2 ) 2 



AM 2 (M 2 + m\-m\- 2Mu x ) 
M(M 2 + m 2 - m 2 2 - 2Muo 1 ) 2 ^ 



+ 00 „ , T 

Mda , 



2vri. 



A 2 A 2 

(- 



2(M 2 - m\ + m\ + 2Mu 2 f 



+ 



16M 6 a; 2 (M 2 - cj 2 + 2Mcj 2 + c^) 2 
(M 2 + m 2 - m\ - 2Mcji) 2 (M 4 + M 2 (mf - m\ - cuf - u\ 
16M 6 uf(M 2 - 2Muj\ + uj\ - ujI) 2 



(ml - m 2 )(3uj 2 - u 2 ) - AMu^m 2 -m 2 + u 2 )) ^ 

16M e ujf(M 2 - 2Mu 1 u 2 ) 2 >' 1 ' 

We have thus evaluated the general expressions for fp (Eq. ff23l ) and Np 
(Eq. fTSHl) ) in the framework of BSE under CIA, with Dirac structure (^75) (7 • 
P)/M introduced in the Hqq vertex function as the leading order structure as 
well as 75, according to our power counting rule. We see that so far the results 
are independent of any model for 4>(q). However, for calculating the numerical 
values of these decay constants one needs to know the constant coefficients Bq 
and Bi which are associated with the Dirac structures 75 and (^75) (7 -P)/M, 
respectively. The relative value, -|Ms a free parameter without any further 
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knowledge of the meson structure in the framework discussed above. As a 
first step, we vary this parameter to see the effect of introducing the Dirac 
covariant ij B j-P/M. We see that at Bi/B = .163, f K = 159.8MeV, which is 
the experimental value of this quantity. For making comparison with results 
of other models and data, we use this value of the ratio Bi/B Q . We also 
vary Bi/B in the range .14 — .17 to demonstrate the dependence on this 
parameter. The results are given in Table I along with those of other models 
and experimental data. It is seen that most of the numerical values of these 
decay constants in BSE under CIA improve when 75(7. P)/M is introduced 
in the vertex function in comparison to the values calculated with only 75. 
Further discussions are in Sec. 4. 

4 Discussions 

In this paper we have first written Hqq vertex function for a pseudo-scalar 
meson employing the power counting rule [I] for the incorporation of various 
Dirac covariants from their complete set. They are incorporated order- by- 
order in powers of inverse of the meson mass. We then calculate fp for 
unequal mass pseudo-scalar mesons (K, D, D$, B) in the framework of 
Bethe-Salpeter Equation under CIA, using the hadron-quark vertex function 
T(q) in Eq. (fT2"j) with the incorporation of the leading order covariants. 
Unequal mass kinematics is employed in this calculation. It is seen that the 
values of Decay constants can improve considerably when the Dirac structure 
275(7 • P)/M is introduced in the vertex function with tuned parameter -|^, 
and come closer to the results of some recent calculations [T2], [131 [191 120] as 
well as agree with experimental results [21 J within error though fp is a bit 
lower (but still within l-a, for details see Table 1). 

In a recent work [12], the Leptonic decay constants fp have been calcu- 
lated for light pseudo-scalar mesons within a ladder-rainbow truncation of 
coupled Dyson-Schwinger and Bethe-Salpeter Equations using Hqq vertex 
function to be a linear combination of four dimensionless orthogonal Dirac 
covariants where each covariant multiplies a scalar amplitude Fi(q 2 , q ■ P, P 2 ) 
for three different parameter sets for effective interactions. The decay con- 
stants for pseudo-scalar meson K calculated in that model [5j [12] for an 
intermediate value of one of the parameter sets: uo = OAGeV, D = 0.93GeV 2 
is fx = 155MeV. In our model with leading covariant7s alone we obtain 
fx = !53.5MeV, which is quite close to this figure. However, here we include 
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the other leading order Dirac covariant 2757 • P/M besides 75, and get fp as 
a function of the parameter ratio Bi/B Q . We then calibrate E>i/B to repro- 
duce the experimental value of Kaon decay constant, fx = 159.8MeV to get 
the value B\/Bq = .163. Using this value of tuned parameter, we then obtain 
f D = 232.78MeV (Expt. 222.6 ± 16.7), f Ds = 295.18MeV(Expt. 294 ± 27), 
which agree with data [2j] within the errors. The decay constant for B- 
meson predicted in our framework is fp = 191.6MeV", which is not far from 
the recent experimental result for fp [22]. Further our model predicts fp s 
value to be around 22% larger than fp value which is roughly consistent with 
the prediction of most of the models which generally predict fp s to be 10% 
- 25% larger than fp (as per recent studies in Ref. [22]). These results are 
demonstrated in Table I. 

Among all the pseudo-scalar mesons, pion enjoys a special status. The 
large difference between the sum of two constituent quark masses and the 
pion mass shows that the quark is far off-shell and q could be the same 
order as the pion mass. So it seems better to incorporate all the four Dirac 
covariants. However, the present experimental condition does not make it 
economical and practical to take into account all four Dirac covariants and 
make a global fitting, for the experimental precision is in different ranks for 
the mesons covered in this paper. On the other hand, from Table I we see 
that the fp value calculated with only the leading order Dirac covariants 75 
and 2757 • P/M rather than all the four Dirac covariants for pion can give the 
results different from data by only about 20%. This indicates the moderate 
contribution from the other two higher order Dirac covariants for pion. Such 
a case can be contrasted with the case of introducing the other LEADING 
order covariant 2757 ■ P/M (besides 75) for three heavy mesons, where we 
see the significant contribution as large as that from 75 term since they are 
equally leading order covariants and the results for f P using these leading 
order covariants for heavier mesons are quite close to their experimental 
values. 

These numerical results for fp obtained in our framework with use of 
leading order covariants 75 and 757 • P/M, as well as those of fy [TJ confirm 
the validity of our power counting rule, according to which the leading order 
covariants should contribute maximum to meson observables, and inspire the 
possibility to investigate the higher order terms in order to get better agree- 
ment between calculations and data (when enough precision is obtained), 
hence help to obtain a better understanding of the hadron inner structure. 
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Table 1: Leptonic decay constants (in MeV) fp in BSE under CIA for 
range of values of ratio Bi/B . The decay constants are calculated from data 
[2T] . The masses of hadrons are also from [5T]. The values of constituent 
quark masses used are set to be m u>d = 300MeV, m s = 540MeK, m c = 
1500MeV, nib = 4500MeK, which is compatible with the other parameters 
fixed from hadron spectrum. Comparisons with results obtained from other 
models are also provided. As discussed in the paper, we calibrate Bi/B to 
Kaon data (the 2nd line), however, if we calibrate Bi/B to pion data the 
results (the 5th line) are as good for the heavier mesons. 
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